function [vN1] = BoundaryDigital(S, K, r, T, tau)

vN1 = exp(-r * (tau)) * normcdf(S, K, 0.01);

